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We present an extension of the multi-purpose Monte-Carlo event generator Sherpa for processes 
in deeply inelastic lepton-nucleon scattering. Hadronic final states in this kinematical setting are 
characterised by the presence of multiple kinematical scales, which were up to now accounted 
for only by specific resummations in individual kinematical regions. An extension of a known 
method for merging truncated parton showers with higher-order tree-level matrix elements allows 
to obtain predictions which are reliable in all kinematical limits. 
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1. Introduction 

Deep-inelastic lepton-nucleon scattering (DIS) offers the possibility to study the structure of 
the nucleon and the dynamics of strong interactions by means of a pointlike probe. It provides 
a clean experimental setting to analyse inclusive quantities like for example the proton structure 
functions. However, the kinematical configurations in DIS are usually very different from those in 
other processes in collider experiments. The virtuality of the exchanged photon tends to be close to 
zero, while final state jets might still have large transverse momenta. Experiments conducted at the 
HERA collider have shown that the large available phase space can lead to a considerable number 
of hard jets, even if the photon virtuality is low. 

This fact poses a certain problem for the simulation of DIS with Monte-Carlo event generators. 
Such programs usually employ parton showers based on the DGLAP equations [1]. It is assumed 
that any scattering process factorises into a core 2 — >■ 2 interaction and a shower evolution, which 
simply "dresses" the hard interaction with softer radiation. As the only hard scale set by the leading- 
order DIS process e q — > e ± q is the photon virtuality, — Q 2 , the probability to produce a jet of 
transverse momentum larger than Q 2 would then vanish. In order to reliably simulate DIS events, 
one must therefore resort to different techniques. We aim at an approach based on combining 
higher-order tree-level matrix elements with the parton shower along the lines of Ref. [2], which is 
consistent with the DGLAP framework. The technical prerequisites for realising this method are 
found in the multi-purpose Monte-Carlo event generator Sherpa [3]. 

2. Event generation technique 

The basic idea of the approach is to separate the phase space into a matrix-element and a 
parton-shower domain through a cut in the phase space of multi-parton processes. The matrix- 
element domain is then supposed to contain hard, well-separated partons only, while the parton- 
shower domain covers the region where resummation effects become important. Throughout the 
hard domain parton-shower emissions are corrected using tree-level matrix elements up to a given 
maximum multiplicity. In the soft domain, the parton shower is applied as is. The separation is 
achieved in terms of a so-called jet criterion, defining the "hardness" and/or the separation of a 
parton with respect to others [2]. Several successful studies of various classes of processes have 
demonstrated the capability of this technique to correctly describe multi-jet final states [4, 5]. 

As pointed out in [2], the above merging algorithm needs to be refined if the scale difference 
between Q 2 and the hardness scale k\ of additional partons is large and negative. In this case, 
logarithmic corrections are not induced by Q 2 /q 2 , but rather by kj/q 2 , where q 2 is the jet resolution 
scale. The production of the virtual photon can then be regarded as an electroweak splitting process, 
attached to a core interaction of type y* j — >■ jj, as depicted in Fig. lb. In the extreme case of very 
hard jets, the core process does not even include the virtual photon, cf. Fig. lc. The main task of 
the algorithm is to correctly identify the "core" interaction underlying a multi-parton process and 
to employ it to define starting conditions for the parton shower evolution [6]. 

Conversely, this idea can be used to lift the restriction on the real-emission phase space at low 
Q 2 . A similar method is in fact employed in Drell-Yan lepton-pair production via y*/Z-exchange, 
where the separation cut Q cut between matrix-element and parton-shower domain is set such that 
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(a) (b) (c) 



Figure 1: Schematic view of three possible core process choices in DIS three-jet production. 
Part (a) corresponds to the most probable core process being the virtual photon exchange, while 
additional hard partons are interpreted as partem shower emissions. Parts (b) and (c) depict config- 
urations, where the most probable core process is the interaction of the virtual photon with a partem 
and a pure QCD 2 — > 2 process, respectively. 



2cut < miii, with miii being the invariant mass of the lepton pair. This choice implies that jets of 
k\ > mui always fall into the matrix element domain. In deep-inelastic-scattering the situation is 
slightly different due to the variable value of Q 2 , which plays the role of miy in the Drell-Yan pair 
production case. The solution can, however, be identical. We choose 

-1/2 

(2-1) 

where Q cut is a fixed value, much like Q cm in the Drell-Yan pair production case and 5dis < 1 
is a constant with lower limit enforced by experimental requirements. Q cut ensures that high-<2 2 , 
medium-is^ B events are described by matrix elements, rather than by the parton shower. At the 
same time, the factor in the square bracket enforces low-<2 2 , high-£y B events to be in the matrix- 
element domain as well, such that the complete available real-emission phase space can be filled. 

3. Comparison with experimental data 

In this section, we show some comparison with experimental data to exemplify the perfor- 
mance of the Monte Carlo simulation. The correct description of the selected measurements is 
quite challenging for the Monte Carlo traditionally used in the analysis of HERA data [7]. 

A crucial observable is given by the inclusive jet cross section, differential with respect to 
Ej B /Q 2 , where Ej.b is the jet transverse energy in the Breit frame. For Ej B /Q 2 > 1 it probes a 
part of the phase space where leading order Monte-Carlo models without the inclusion of low-.* 
effects are bound to fail in their description of jet spectra. Figure 2a shows that the Monte-Carlo 
prediction gradually improves with a growing number of final-state partons in the hard matrix 
elements. The uncertainties associated with a variation of the intrinsic parameters of the merging 
algorithm are shown in Fig. 2b. 

It is interesting to investigate jet properties in some more detail. Figure 3 displays rapidity 
spectra of the forward jet in di-jet production for various regions of Q 2 . As for the case of inclusive 
jet production we observe a good description of the corresponding HI data. 
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Figure 2: The inclusive jet cross section as a function of Ej B /Q 2 in bins of Tf]i a b, measured by the 
HI Collaboration [8]. Ej B is the jet transverse energy in the Breit frame, while T}/^ denotes the jet 
rapidity in the laboratory frame. Part (a) displays the influence of the maximum partem multiplicity, 
A^ max , from hard matrix elements. We show the uncertainty originating from varying 5dis between 
0.5 and 0.7 (light grey band) and from varying Q cut between 3 GeV and 9 GeV (dark grey band) in 
part (b). 



4. Conclusions 

The Sherpa event-generation framework has been extended to describe hadronic final states 
in deep-inelastic lepton-nucleon scattering processes. The simulation is based on merging higher- 
order tree-level matrix elements with a partem shower. When applying this technique to DIS pro- 
cesses, it is vital to correctly identify the core interaction, which can be either electron-quark scat- 
tering, photon-parton scattering or a partonic 2 — > 2 interaction, depending on the final state kine- 
matics. The particular kinematical situation in DIS also requires to choose appropriate merging 
scales, depending on the photon virtuality Q 2 . By doing so, we obtain a reliable description of DIS 
in all kinematical regions, including for example high-is^, low-<2 2 processes. 
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Figure 3: The di-jet cross section as a function of T7f w ,iab in bins of Q 2 , measured by the HI 
Collaboration [9]. 
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